Chemical freeze-out in relativistic ^4+^4 collisions: 
is it close to the QGP? 
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Abstract 



Preliminary experimental data for particle number ratios in the collisions of 



^ i Au+Au at the BNL AGS (11 A GeV/c) and Pb+Pb at the CERN SPS (160,4 

GeV/c) are analyzed in a thermodynamically consistent hadron gas model 
with excluded volume. Large values of temperature, T = 140-185 MeV, and 
baryonic chemical potential, = 590-270 MeV, close to the boundary of 



the quark-gluon plasma phase are found from fitting the data. This seems to 



indicate that the energy density at the chemical freezeout is tremendous which 

would be indeed the case for the point-like hadrons. However, a self-consistent 

treatment of the van der Waals excluded volume reveals much smaller energy 

densities which are very far below a lowest limit estimate of the quark-gluon 

plasma energy density. 
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The quantitative results of quantum chromo dynamics for the equation of state of 
strongly interacting matter have been obtained mostly in the high energy density region, 
e > 2 GeV/fm 3 , where one expects a weakly interacting gas of quarks and gluons — the 
quark-gluon plasma (QGP). The region of small e corresponds to the hadron gas (HG) 
phase and can only be treated today in terms of the hadron degrees of freedom. Assuming 
a local thermodynamical equilibrium at the final freeze-out stage of nucleus-nucleus (A+A) 
collisions, one can estimate the particle number ratios without detailed knowledge of the 
complicated system evolution. Hadron abundances and ratios have been suggested as pos- 
sible signatures for exotic states and phase transitions in dense nuclear matter (see, e.g., 
Refs. [1-5]). 

Preliminary data for A+A collisions with truly heavy beams have recently become avail- 
able: Au+Au at 1L4 GeV/c at the BNL AGS and Pb+Pb at 160A GeV/c at the CERN 
SPS 0. A systematic analysis of these data could yield clues to whether a short-lived QGP 
phase with high e exists during the hot and dense stage of these reactions. Recent fits of 
particle number ratios [7-12] at AGS and SPS energies in the framework of "thermal model" 
have led to large values of temperature, T, and baryonic chemical potential, fib, close to the 
boundary of the QGP phase |TT| , |T2"|] . The aim of the present paper is two-fold. First, we 



clarify the notion of the "thermal model" used in the literature. Second, we present our 
results for Au+Au (AGS) and Pb+Pb (SPS) collisions and examine whether the chemical 
freeze-out states are really close to the QGP. We find large values of the chemical freeze-out 
parameters T = 140-185 MeV and fib = 590-270 MeV. The values of the energy density 
e and baryonic number density are, however, strongly dependent on the specific ther- 
mal model formulation. In the ideal HG, the chemical freeze-out values of e and rib for SPS 
A+A collisions are indeed close to the QGP estimates. These values become, however, much 
smaller in thermodynamically consistent HG model with excluded volume. 

The name "thermal model" has always been used in the literature when calculations of 
particle number ratios in A+A collisions are done with parameters T and fib- We stress, 
however, that different people used, in fact, very different "thermal models". It seems 



natural to use the ideal HG as a thermal model at the freeze-out stage. However, such an 
ideal gas model becomes inadequate in high-energy A+A collisions. The chemical freeze- 
out parameters T and fib obtained from fitting the particle number ratios at AGS and 
especially SPS energies lead to artificially high particle number densities which contradict 
the assumed picture of non-interacting hadrons. The van der Waals ( VDW) excluded volume 
procedure has recently been used in a number of thermal model calculations. We follow the 
excluded volume procedure of Ref. |R|. The "thermal models" of Refs. [7-12] also include 
VDW "corrections", but in ad hoc thermodynamically inconsistent ways. Besides, in the 
formulation of Refs. [7-10], particle ratios always remain the same as in the ideal gas for 
any choice of proper hadron volumes. 

For a fixed particle number N, the VDW excluded volume procedure amounts to the 
substitution of the volume V by V — vN, where v is the parameter corresponding to the 
proper volume of the particle. Note that this VDW procedure, interpreted in statistical 
mechanics as an approximation to the gas of hard-sphere particles with radius r, requires 
that the volume parameter v is equal to the "hard-core particle volume" , tvrr 3 , multiplied by 



a factor of 4 [TJj]. An extension of the excluded volume procedure for several particle species 
i = 1, h, has been done with the bf ansatz of the substitution V — > V — Ya=i v i^i . The 
pressure function in the grand canonical formulation is given then by the equation (13| : 
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with di the number of internal degrees of freedom (degeneracy), //j the particle chemical 
potential and the mass. Particle number density for the i-th species is calculated from 
Eq. © as 



with 



Hi = Hi - Vi p{T, hi, fj, h ) , i = l,...,h. (3) 



The energy density of the HG with VDW repulsion is 



(4) 



1 + Ej=i vj nf(T,fij) ' 
where n- d , e* d are the ideal gas quantities. The VDW repulsion leads to strong suppression 
of particle number densities and therefore all thermodynamical functions in the HG with 
the VDW repulsion become much smaller than those in the ideal gas at the same T and 
Hi, ...,Hh (see Ref. [|nj for details). 



Particle chemical potentials Hi regulate the values of conserved charges. For simplicity 
we neglect the effects of non-zero electrical chemical potential which were considered in 



Ref. [ I6|,r7| . The chemical potential of the i-th particle can then be written as 

Hi = hHb + SiHs , (5) 

in terms of baryonic chemical potential Hb an d strange chemical potential Hsi where 6, and 
Si are the corresponding baryonic number and strangeness of the i-th particle. The hadronic 
gas state is defined by two independent thermodynamical parameters, T and Hb- The strange 
chemical potential Hs{T,Hb) is determined from the requirement of zero strangeness 

h 

n s (T,Hb,Hs) = ^2sini(T,Hi) = 0. (6) 
i=i 

For the SPS A+A data it has been observed in Refs. |8||| that their "thermal model" 
(for particle number ratios it is, in fact, equivalent to the ideal gas model) is unable, within 
a single set of freeze-out parameters, to reproduce simultaneously the strange particles and 
anti-baryon to baryon ratios together with ratios where pions are involved. Experimental 
pion to nucleon ratio and ratios of pions to other hadrons are larger than the ideal gas 
predictions. Several mechanisms, including the chemical non-equilibrium with Hn > in the 
pion subsystem and the possibility of straight QGP hadronization P,p|,[l8| were proposed, but 



no satisfactory solution has been found. Another possibility with an incomplete strangeness 
chemical equilibrium and parameter 7 S < 1 has been recently studied in Ref . [HJ . We follow 
the HG model with complete chemical equilibrium. To enlarge the number of pions in our 
VDW HG we assume that the pion proper volume parameter tv is smaller than the proper 
volume parameters for other hadrons which we put to be equal to v. In the case where not 
all of the Vi's are equal, the hadron volume parameters V{ do influence the particle number 
ratios through the modification of Eq. (g) in the particle chemical potentials, required by 
thermodynamical self-consistency. The effect is quite evident: hadrons which take up less 
space (i.e., smaller values oivi) are preferably produced. The particle number ratios of those 
small hadrons to larger ones increase as compared with the ideal gas results. We assume 
different hard-core radii: r n for pions and r for all other hadrons (r > rv), with Vi = 4- ^ rf. 

Baryon and meson resonances are of great importance for the measured particle number 
ratios at AGS and SPS energies. All known resonance states with masses up to 2 GeV are 
included in our calculations with subsequent resonance decays to stable hadrons. We use 
the compilation of experimental data for particle number ratios in Au+Au at the BNL AGS 
(11A GeV/c) and Pb+Pb at the CERN SPS (160 A GeV/c) which were presented by Stachel 
at QM'96 (see Ref. \T2\ and references therein). Our fits within the VDW HG model is shown 



in Fig. 1 (see also |y|). From fitting the data we found the following model parameters: 

AGS : T = 140 MeV , fi b ^ 590 MeV , (7) 



SPS : T = 185 MeV , fi b = 270 MeV . (8) 

There are a number sets of possible values for hard-core radii which give the same hadron 
ratios in Fig. 1. If we put r n = 0, then one finds r = 0.50 fm and r = 0.46 fm for AGS and 
SPS cases, respectively. On the other hand, the AGS and SPS data in Fig. 1 can be fitted 
simultaneously with following set of 

rv = 0.62 f m , r = 0.8 fm . (9) 
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Due to TV < r the ratios of direct pions to other hadrons increase by a factor of 2.0 for AGS 
(0) and 2.6 for SPS (^) in comparison to the case tv = r (in this later case all ratios are 
very close to their ideal gas values). It results in the increase of total pion to other hadron 
ratios about 25% for AGS © and 33% for SPS (§. 

The calculated energy density e for the above T, fi b values are 

AGS : e id 0.72 GeV/fm 3 , e VDW = 0.25 - 0.08 GeV/fm 3 , (10) 

SPS: e id 1.58 GeV/fm 3 , e yw = 0.43 - 0.10 GeV/fm 3 , (11) 

for the ideal HG (i.e. tv = r = 0) and VDW HG, respectively. The corresponding values 
for the baryonic number density n& are 

AGS: nl d = 0.40 fin" 3 , n\ DW = 0.13 - 0.04 far 3 , (12) 

SPS: nl d = 0.40 far 3 , n\ DW = 0.10 - 0.04 far 3 . (13) 

The higher values of e VDW and n( DW in the above equations are obtained for tv = 0, 
r = 0.50 fm (AGS) and r n = 0, r = 0.46 fm (SPS) while lower ones are found with 
parameters of Eq. (|9|) which fit simultaneously the AGS and SPS data. Note also that 
Eq. (3) gives an additional suppression by a factor of exp(— v,- t p/T) for all particle number 
and energy densities of i-th particle in comparison to the inconsistent VDW gas treatment 
[7-10]. In the region of different T, tv and r considered above it leads to a suppression 
factor between 2 and 6 for the total energy density and baryonic number density. 

To compare the energy densities of HG (0,0) with those of the QGP we use equation 
of state of the bag model, i.e., an ideal gas of u,d,s quarks, antiquarks and gluons with 
nonperturbative "vacuum pressure", B = 400 MeV fm -3 (see, e.g. [0]). The chemical po- 
tentials of quarks and antiquarks are given by the general equation (^) with hi = 1/3 (—1/3) 
for all quarks (antiquarks), Sj = for u, d quarks and antiquarks, Sj = — 1 (+1) for s quark 
(antiquark). The consequence of zero strangeness requirement © is quite simple in the 
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QGP. It leads to equal values for the chemical potentials of s and s and hence both should 
be equal to zero. Thermodynamical functions within zero quark mass approximation can 
be easily calculated: 

n? = 5 r2 + 8^ (16) 

The lowest estimate of the energy density for the QGP can be obtained by putting 

p Q (T^ b ) = 0. (17) 

which defines a curve in the /i&-T plane as shown in Fig. 2. The values of (pl|) (s® = 4B) 
and nf ( |i~6|) along this curve are shown in Fig. 3. Note that the curve in the plane 
in Fig. 2 is not realistic for the HG-QGP phase transition. This is just a lower estimate 
for the fjbb-T boundary of such a transition. The values of T and given in Eqs. (|7],[8|) are 
also shown as points in Fig. 2: AGS (0) square-point is very close to the QGP boundary 
and SPS (|8|) circle-point is in the "QGP phase" . Our Fig. 2 is similar to the corresponding 



figures in Refs. ||lT| , p!2| . To compare the energy densities of HG and QGP states one should 



remember that the energy density values are very different in the ideal HG and VDW HG. 



The ideal HG model was used recently to describe particle number ratios in AGS |17j and 



SPS |L9| A+A collisions. Note that large values of T = 193 MeV and fj. b = 234 MeV, close 
to our estimates @, were found in Ref. Jl9| for Pb+Pb (SPS). Fig. 3 shows that ideal HG 
value of e td for the SPS parameters (|8[) is indeed rather close to the QGP phase. However, 
the VDW HG values e VDW for the same T and [i b are very much smaller than e %d . Ideal 
HG and VDW HG models lead to very different pictures of the chemical freezeout for A+A 
collisions at SPS energies: e id = e Q , but e VDW « e Q and n\ DW « rtf . We stress that 
very large values of e ld and n td are in an inevitable logic contradiction with the assumed 
picture of non-interacting hadrons. 



In summary, we have analyzed the preliminary data for particle number ratios measured 
in the collisions of Au+Au at the BNL AGS (11 A GeV/c) and Pb+Pb at the CERN SPS 
(16CL4 GeV/c) with a thermodynamically consistent hadron gas model with excluded vol- 
ume. Even though large values of temperature, T = 140-185 MeV, and baryonic chemical 
potential, [i h = 590-270 MeV, close to the boundary of the quark-gluon plasma phase are 
found from fitting the data, the energy densities obtained always lie far below a lowest limit 
estimate of the QGP energy density because of the VDW suppression effects. Such a picture 
is consistent with a large energy density change during a HG-QGP phase transition. 
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FIGURES 



Fig. 1: Points are the preliminary experimental data for the particle number ratios (see 



Ref. [I2| and references therein) for Au+Au AGS and Pb+Pb SPS collisions (in the lower 
and upper part of the figure respectively). The short horizontal lines are the model fit with 
T = 140 MeV, fi b ^ 590 MeV (AGS), and T = 185 MeV, fi b S 270 MeV (SPS). In both 
cases 7v = 0.62 fm, r = 0.8 fm. 

Fig. 2: The solid line is the curve p®(T, //&) = fli~7D in the //&-T plane. The square and 
circle are the chemical freeze-out points of (0) and @, respectively. 

Fig. 3: The straight solid line gives the lowest energy density estimate, = 4B = 
1.6 GeV/fm 3 , for the QGP. It corresponds to the curve p®(T, //&) = (|1~7|) of Fig. 2 in the 
n b -e plane. The largest value in the x-axis, n b = 1.08 fm~ 3 , corresponds to the point T = 
of the curve in Fig. 2. The open square and circle show the ideal HG values of e id , n l b d for 
AGS (0),© and SPS 0,(0). The solid lines show the sets of the VDW HG values 
£ vdw ^ n VDW f or AQg anc [ sps. The solid square and circle are the VDW HG values at 
AGS (0) and SPS (g) with hard-core radii (§. 
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